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We examine the effects of hole— phonon interaction on the formation of hole pairs in the 
2D Holstein t-J model. Using finite— lattice diagonalization techniques, we present exact 
results for the two— hole binding energy and the s- and d-wave pairing susceptibilities. 

The interplay of electronic and lattice degrees of freedom in strongly electron- 
correlated systems is now attracting a lot of attention. This interest is partially due 
to the prominent role of the electron-phonon (EP) interaction in several transition 
metal oxides with strong Coulomb correlations, such as the colossal magnetore- 
sistive manganites or the charge-ordered (insulating) nickelates. For the high-Tc 
superconducting cuprates, very recent experiments demonstrate the relevance of the 
coupling between the charge carriers and the lattice dynamics as wellS. 

For the sake of simplicity, we describe the basic interactions in the latter_aystems 
by an effective single-band Hamiltonian, the so-called Holstein t-J modella 

n - H,_j - V^oY. (^1 + + hoJo {bib, + 1) , (1) 

i i 

which contains besides nearest-neighbour hole transfer (t) and antiferromagnetic 
(AFM) spin exchange (J) on a square lattice, the couphng of doped holes (hi) to a 
dispersionsless optical phonon mode (representing, e.g., local apical-oxygen breath- 
ing vibrations). Here Sp is the hole-phonon coupling strength and hujo denotes 
the phonon frequencv.. The single-particle excitations of the model (1) have been 
studied numerically B'eI; the main result is that in the presence of strong AFM spin 
correlations even a moderate EP coupling can cause polaronic effects. 

In this contribution, we focus on the two-hole subspace in order to comment on 
hole-pair formation. Employing the Lanczos algoritkp in combination with a well- 
controlled truncation of the phononic Hilbert spaceHiJ, we are able to calculate both 
ground-state and dynamical properties preserving the full dynamics and quantum 
nature of phonons. It is worth noticing, that our multi-mode treatment of the 
phonons differs significantly from the one-phonon calculation performed in Ref. 6 
for the t-J model coupled to oxygen breathing or buckling modes. 

The dynamical pair spectral function can be written as 

A,,{u;) = ^ |(f I I 5 - (Si^-^) - , (2) 

n 

where Aa = -^J2i;S=±x,±y^a{5)ciiCt+si, with Tai±x) = 1 and Ta{±y) = -1 
[+1] for d-wave [s-wave] pairing. 

Numerical results for the spectral functions of the pairing operators Aa are 
shown in Fig. 1 for exchange interactions J = 0.4 (left column) and J = 0.1 (right 
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Figure 1: Dynamical pair spectral function A2h{'^) (dotted lines) and integrated spectral weight 
J^2h{'^) = duj' A2h{^') (bold lines) calculated for the 2D Holstein t-J model on a ten- 

site square lattice with periodic boundary conditions at Huiq = 0.8. Depending on the model 
parameters J and Sp, results arc presented only for those hole-pair wave functions Izia!^','^') 
which have a finite overlap with |!^'g^ ^''). Note that for the the pure t— J model the symmetry 
of the two— hole (K = 0) ground state is changed from d— wave {B/C4) to s-wave (A/C4) at 
Jc = 0.2001 {N = 10). On the other hand, the (symmetry) change in the spectra, observed by 
comparing the plots for £p = 0.5 and 3.5 at J = 0.1, is driven by the EP coupling (£p,c — 3.3). 
All energies are measured in units of t. 
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column) corresponding to two different regimes in the pure t-J model (see upper 
panels). For J — 0.4, the d-wave pair spectrum exhibits a well-separated low- 
energy peak containing an appreciable amount of spectral weight which grows if 
J is enhanced (cf. the inset of Fig. 2). Since the rest of the spectrum becomes 
incoherent with increasing lattice size iV, the dominant peak at the bottom of 
the spectrum has been taken as signature of a d-wave quasiparticle bound state Q. 
By contrast, the s-wave spectrum shows no such quasiparticle-like excitation. In 
the weak EP coupling case, the main features of the = spectra are preserved, 
although, of course, additional phonon satellite structures appear (cf. the discussion 
of Fig. 3 below). The situation is drastically different in the strong-coupling regime. 
Here a strong mixing of electron and phonon degrees of freedom takes place and less 
mobile (bi)polaronic charge carriers emerge. The polaronic self-trapping transition 
is accompanied by a dramatic reduction of the coherent band width and, as a result, 
the AFM spin interaction becomes much more effective. Therefore, at e-p — 3.5, the 
spectrum for J = 0.1 looks very similar to that for J = 0.4. 

The relative spectral weight, Z2h, located in the lowest pole of ^2/1, is plotted 
in Fig. 2 (left panel). Obviously, we observe a strong suppression of the d-wave 
quasiparticle residue with increasing EP coupling. However, Z^h gives only a mea- 
sure of the "electrjonic" contribution to the d-wave bound state. In fact, according 
to previous worku, composite pair operators A^^ properly dressed by a phonon 
cloud, give large quasiparticle weights in A2h- That the Holstein EP interaction 
may stabilize a bound state of two holes is clearly demonstrated by the behaviour 
of the binding energy i?^, which has been calculated for the larger 4x4 lattice in 
order to reduce finite-size effects (see the right panel of Fig. 2). Whereas the hole 
attraction {E\ < 0) is less affected in the anti-adiabatic limit, the EP coupling 
promotes the pairing correlations between two holes in the adiabatic regime due to 
subtle retardation effects. This is to be contrasted with the findings for a coupling 
of the holes to the in-plane oxygen breathing mode which leads to an hole repulsion 
for all frequencies and EP interaction strengths^. 




Figure 2: Two-hole "quasiparticle weight" Z2h = K'^'o^ I I '^'q^')! VK^^'o^' l^l^al^'o'^^ )\ 
(left panel; A'' = 10) and hole "binding energy" Eg = iJg^"^' + -Bq'^'' — 2Eg'^~^^ (right panel; 
A'^ = 16 ) shown as a function of EP coupling strength Ep at various phonon frequencies hajQ. 
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Let us discuss the weak-coupling case in more detail. Fig. 3 presents the spec- 
tral decomposition of the d-wave pairing operator at huJo — 0.1 and 3.0. In the 
high phonon frequency (anti-adiabatic) regime, the pairing susceptibility behaves 
qualitatively in a similar way as the £p = limit, in particular the low-energy part 
of the spectrum is given by purely electronic resonances. If the phonon frequency is 
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Figure 3: Left panel: d— wave pair spectra A2h £p = 0-5 a^nd huj = 0.1 (inset: hio = 3.0). Right 
panel: integrated spectral weight M2h{^) in the low— energy part of A2h for various ep (inset: 
phonon— weight function |c™p in the two-hole ground state for the same parameters) 



much smaller than the electronic gaps, we found series of predominantly phononic 
side bands, being separated by Tlwq and roughly centered around the positions of 
the electronic excitations. The relative weights of these 5-like peaks can be de- 
duced from the corresponding jumps in A/2/1 (w) depicted for different Sp in the right 
panel of Fig. 3. Focusing on the lowest set of phonon sub-bands, it is interesting 
to note, that the weight of the zero-phonon state (|c°p) in the first excited states 
(l^l'n^ ^''), n ~ 1, 2, . . .), which is measured by J\f2hi'-^), is approximately the same 
as the weight of the m-phonon states (|cq'|'^) in the ground state {\^^'^^)). The 
definition of the coefficients |cq'P is given in Refs. 4,5. A qualitative understanding 
of this fact can be obtained from the study of the independent boson modelB, where 
we can show exactly that |c™p = \c^\'^ holdsu. 

To summarize, our exact diagonalization studies of the 2D Holstein t-J model 
give evidence for significant EP coupling effects. Most notably, the hole pairing in 
the t-J model may be stabilized by a dynamical (Holstein) hole-phonon interaction. 
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